Initial state dependence of the quench dynamics in integrable quantum systems. 

II. Thermal states 
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We study properties of isolated integrable quantum systems after a sudden quench starting from 
thermal states. We show that, even if the system is initially in thermal equilibrium at finite tem- 
perature, the diagonal entropy after a quench remains a fraction of the entropy in the generalized 
ensembles introduced to describe integrable systems after relaxation. The latter is also, in general, 
different from the entropy in thermal equilibrium. Furthermore, we examine the difference between 
the distribution of conserved quantities in the thermal and generalized ensembles after a quench and 
show that they are also, in general, different from each other. This explains why these systems fail 
to thermalize. A finite size scaling analysis is presented for each quantity, which allows us making 
predictions for thermodynamically large lattice sizes. 
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I. INTRODUCTION 

The recent renewed interest in the theoretical study of 
the relaxation dynamics and description after relaxation 
of isolated many-body quantum systems is largely mo- 
tivated by extraordinary advances in experiments with 
ultracold atomic gases, which provide highly tunable set 
ups with very high degree of isolation from the environ- 
ment [241 . Among many remarkable experimental re- 
sults, one can mention the collapse and revival of matter- 
wave coherence after the sudden quench from a super- 
fluid to a Mott insulator [l|, the finding that one- 
dimensional (ID) Bose gases can relax to nonthermal 
equilibrium states [H, and the almost perfect agreement 
between the nonequilibrium dynamics observed in experi- 
ments with ID lattice bosons and numerically exact simu- 
lations of the unitary dynamics using the time-dependent 
renormalization group approach (sj . 

Motivated by the findings in Ref. Q , intensive theoret- 
ical efforts have been devoted to understanding the de- 
scription after relaxation of isolated integrable systems 
following a sudden change in a Hamiltonian parameter 
(sudden quench). In quenches between integrable sys- 
tems, in which the initial state is an eigenstate of an 
integrable Hamiltonian (usually the ground state), many 
studies have shown that observables do not relax to the 
thermal predictions [6l-l23| . Instead, they relax to the pre- 
dictions of generalized Gibbs ensembles (GGEs), which 
take into account the presence of nontrivial sets of con- 
served quantities j6l-l25|. The GGE maximizes the en- 
tropy while taking into account the constraints imposed 
by those conserved quantities @, [H, [13] ■ 

Microscopically, the GGE has recently been shown to 
work for ID lattice hard-core bosons [1^ because of a 
generalization of the eigenstate thermalization hypothe- 
sis H, [1^ . The idea is that the eigenstates of the 
final Hamiltonian that overlap with the initial state have 
very similar expectation values of the conserved quanti- 
ties (properly accounted for in the GGE) and very similar 
expectation values of few-body observables, i.e., ETH is 



valid within that restricted set of eigenstates [l9| . Fur- 
thermore, for quenches very close to integrable points, 
there is a separation of time scales in which the system 
usually relaxes to a prethermalized state before thermal- 
izing [30l433| . Observables in that prethermalized state 
have been shown to be described by a GGE [13] ■ Inter- 
estingly, for quenches from eigenstates of nonintegrable 
systems to an integrable point, it has been argued that 
thermalization can occur because the initial state pro- 
vides an unbiased sampling of the eigenstates in the final 
Hamiltonian, very much as the thermal ensembles do [ssj . 

As mentioned before, most works that have directly 
probed the nonequilibrium dynamics of isolated inte- 
grable systems have focused on initial states that are 
eigenstates (usually the ground state) of the Hamilto- 
nian before the quench. This opens questions as to how 
general the conclusions obtained in those studies are for 
more generic finite-temperature systems. Thermal initial 
states were considered for quenches in the quantum Ising 
model in Ref. [s^. There, it was shown that for suffi- 
ciently high initial temperatures nearly thermal distribu- 
tions occur for the conserved quantities after a sudden 
quench. Since away from the critical point the quantum 
Ising model is gapped, "sufficiently high" in this case 
means higher than the values of the relevant gaps. 

In this work, we study what happens for quenches 
starting from equilibrium finite-temperature states in ID 
lattice hard-core bosons (the XY model). Within this 
model, a local chemical potential (a site-dependent z 
magnetic field) can be used to generate gaps in the spec- 
trum. We focus on the resulting entropy and distribu- 
tion of conserved quantities after a sudden quench and 
compare them to the predictions of thermal ensembles 
and the GGE relevant to the final system in equilibrium. 
Understanding the outcome of the relaxation dynamics 
for initial states at finite temperature is particularly im- 
portant to address current ultracold gases experiments, 
where quenches are performed starting with the gas at 
some effective finite temperature. 

In the previous and closely related work, we studied 
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the same quantities and quenches similar to the ones 
considered here but starting from the ground state of 
the initial Hamiltonian [2l|. We showed that, if the ini- 
tial state was the ground state of a half-filled system in 
a period 2 superlattice (an insulating state), the distri- 
bution of conserved quantities after the superlattice was 
turned off approached that in thermal equilibrium upon 
increasing the superlattice strength. (An understanding 
of this phenomenon in terms of bipartite entanglement 
was recently provided in Ref. [l^l-) At the same time, 
the entropy of the GGE approached that of the ther- 
mal ensemble (grand-canonical ensemble; GE). However, 
contrary to what happens in nonintegrable systems, the 
difference between the GGE and the thermal ensemble 
predictions for the entropy and the conserved quantities 
did not vanish in the thermodynamic limit for any finite 
strength of the superlattice. For all other quenches con- 
sidered in that work, the predictions of the GGE and the 
thermal ensembles were quantitatively and qualitatively 
different. 

Here we show that following a sudden quench and 
upon increasing the temperature of the initial state, sim- 
ilarly to what was found in Ref. [2l| when increasing the 
strength of the initial superlattice potential and consis- 
tent with the findings in Ref. (36| for the quantum Ising 
model, the distribution of conserved quantities and the 
entropy in the GGE approach the thermal predictions. 
However, this should not be confused with thermaliza- 
tion, as, for any given initial finite temperature and finite 
Hamiltonian parameters before and after the quench, the 
differences between the two ensembles remain finite in 
the thermodynamic limit. Hence, our results show that 
for quenches between integrable systems, there is no fun- 
damental difference between starting with an eigenstate 
of the Hamiltonian and starting with a state in thermal 
equilibrium. Thermalization does not generally occur in 
either case. 

The exposition is organized as follows. In Sec. [Hi we 
introduce the model of interest, and the definition of 
the statistical ensembles and the observables studied. In 
Sec. IIII[ we discuss the behavior of the weights of the 
eigenstates of the final Hamiltonian in various ensem- 
bles. The energy density distributions are also studied 
in that section. Subsequently, in Sec. IIVI we study var- 
ious entropies and perform a series of scaling analysis 
to assess finite-size effects as well as the influence of the 
control parameters in the properties of the system after 
the quench. A similar study, but for the behavior of the 
conserved quantities, is presented in Sec. [V] Finally, our 
results are summarized in Sec.lVIl 



integrable model with Hamiltonian 



II. MODEL, ENSEMBLES, AND OBSERVABLES 

We focus on the nonequilibrium properties of ID lattice 
hard-core bosons following a sudden quench. This is an 
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H = -tJ2(^lk+i + H.c.) + A^(-l) 



(1) 



where t is the hopping parameter and A is the strength 
of a local alternating (superlattice) potential. We con- 
sider lattices with L sites and open boundary condi- 
tions. The hard-core boson creation (annihilation) op- 
erators are denoted bl(bi), and the number operator 
fii = blbi. In addition to the bosonic commutation re- 
lations [61,6]] = Sij, there applies a constraint that sup- 
presses multiply occupancy of the lattice sites in all phys- 
ical states, i.e., b]^ = fef = 0. We note that, in the fol- 
lowing, h = 1, ks = 1, and the hopping energy is set to 
t = 1 (our unit of energy) . 

To be solved exactly, this model can first be mapped 
onto the spin- 1/2 XY Hamiltonian through the Holstein- 
Primakoff transformation [38[ and, subsequently, via the 
Jordan- Wigner transformation [39| . to noninteracting 
fermions. Instead of diagonalizing the full many-body 
Hamiltonian, one can then write each many-body eigen- 
state in the form of a fermionic Slater determinant with 
appropriate modifications following the mapping rules. 
These Slater determinants are constructed as products of 
single-particle eigenstates of the noninteracting fermionic 
Hamiltonian. Utilizing properties of Slater determinants, 
one can com put e exactly all one-particle [13, El and two- 
particle (42!. |43|| observables in the eigenstates, at finite 
temperature in the GE , and out of equilibrium [isj . 

Before the quench, our system is assumed to have N 
particles and to be in contact with a thermal reservoir 
at finite temperature T in equilibrium. At time r = 0, 
we disconnect the system from the reservoir and change 
the strength of the alternating potential from its initial 
value Aj to its final value Ap {Hi — > Hp)- The initial 
thermal state, which is a mixed state, can be expressed 
in the basis of the many-body eigenstates j^*^) of the 
initial Hamiltonian Hj, Hi\^i) = E^^fi). The many- 
body density matrix can be written as 
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1 



(2) 



where Zr = e ^^/T ^ 'pj-^g time evolution of the den- 
sity matrix is given by 



/5(r) = ^E^"'^"^''l<W)(*iMI' 



(3) 



in which the time-evolving many-body eigenstates of 
the initial Hamiltonian Y^a{T)) can be written in terms 
of the many-body eigenstates of the final Hamiltonian, 



E 



(4) 
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Following Eq. (|4]) , the infinite time average of Eq. ([3|) can 
be written in the form of a diagonal density matrix, which 
corresponds to the so-called diagonal ensemble (DE) Q , 

~p= \\ui -( dTp{T)=pT,^=Y,Wp\^){m 

t'->oo r /n ^ — ' 

(5) 

where we have assumed that degeneracies, if present, are 
irrelevant, and have defined 



Here, Wfi corresponds to the weight of state 
the DE. As we show below, is strongly dependent 
of the initial state temperature and the quench proto- 
col. In Eq. the overlaps between eigenstates of the 
initial and those of the final Hamiltonian (^'^I^'q) are 
evaluated numerically as the determinant of the product 
of two matrices, with each matrix composed of the ma- 
trix elements of the Slater determinant representing the 
eigenstates 2l|. Note that the dimension of the Hilbert 
space for our model is d = (^y) , and the computation 
time for the entire set of DE weights is proportional to 
cP ^ i.e., it is exponential in the system size. 

The expectation value of any observable O after relax- 
ation, if relaxation takes place, is equal to the infinite- 
time average O and can be calculated as Tr[/5DE O] [loj . 
Following the definition of the von Neumann entropy 
S = — Tr[/51n/5], one can also calculate the entropy of 
the DE as 



5, 



DE 



Y^WpHWp). 



(7) 



This entropy has recently been shown to satisfy the re- 
quired properties of a thermodynamic entropy, namely, 
it increases when a system in equilibrium is taken out of 
equilibrium, it satisfies the fundamental thermodynamic 
relation [4^, and it is additive and equal (up to subex- 
tensive corrections) to the entropy of thermal ensembles 
used to describe generic systems after relaxation [47| . 

Here, the DE is compared with various statistical 
ensembles. If the number of particles N in the sys- 
tem is kept fixed and the energy is allowed to fluctu- 
ate about a mean value determined by the initial state 
Ej = Tt[piHp\, the relevant ensemble to compare with 
is the canonical ensemble (CE). The density matrix in 
this case has the form 



1 

^CE 



(8) 



where Zce = X^^j ^ ^'^ ^^^^ ^^-^ Tqe is taken so that 
Tr[/5cE^_F] = El. Similarly, for a system in which not 
only the energy but also the number of particles is allowed 
to fluctuate, the relevant ensemble is the GE, for which 



the density matrix has the form 
1 



Pge — 



^GE 



^ ^-{E'f< -mA^I I^F^^^ 
/3 



(9) 



where Zge = Z^a^ 



13 )/TgE g^^^ J\jF gjygj^ 



Nq\^^) {N is the operator for the total num- 
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ber of particles). Note that, in Eq. the sum runs over 
the entire set of many-body eigenstates of Hp with all 
possible particle numbers. Tqe and need to be chosen 
so that Tr[pGE^F] = Ej and Tr[pGE-/V] = N. The en- 
tropies corresponding to each of those ensembles can be 
directly calculated using von Neumann's definition and 
yield 



iScE = In ZcE - 
'S'ge = In Zqe 



Ei 

TcE ' 

El- pN 

?GE 



(10) 

(11) 



These two entropies agree with each other in the thermo- 
dynamic limit, for any given value of N and E^ and the 
same is true for the predictions of both ensembles (and 
the microcanonical ensemble) for additive observables. 
In the GE, because of the mapping between hard-core 
bosons and fermions, Zqe can be calculated in terms of 
the single-particle eigenenergies of the final Hamiltonian 
e„, and Tge and p [4J]: 



^GE = 



n[ 



(12) 



In isolated integrable systems, however, the presence 
of nontrivial sets of constants of motion constrain the 
dynamics and lead to expectation values of observables 
after relaxation that arc different from the predictions of 
thermal ensembles. The GGE is then the natural choice 
to describe such systems 0. The many-body density 
matrix of the GGE can be written as 



1 



Pgge 



GGE 



(13) 



where ^gge = Tr[e and, for hard-core bosons, 

{/„} are nothing but the projection operators to the 
single-particle eigenstates of Hp- (The fact that these 
quantities are conserved is evident if we write the final 
Hamiltonian as Hp = £nln-) {^n} are Lagrange mul- 
tipliers, which need to be selected to meet the initial 
conditions Tr[/5//„] = Tr[/3QGE^ji]- For hard-core bosons, 
they can be computed as |6| 



In 



1 - Tr[p,/„] 



Tr[pil„] 

The GGE entropy can then be written as 

Sgge = In ^GGE + ^ A„Tr[p//„] 



(14) 



(15) 
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In this work, we focus on two types of quenches; the 
first one is the turning-off of the superlattice potential 
Aj ^ 0, Ap — 0, and the second one is the reverse of 
the first one, namely, the turning-on of the superlattice 
potential Aj = 0, Ap ^ 0. In all cases, we start from a 
finite-temperature mixed state. 



III. WEIGHTS IN THE ENSEMBLES 

We are first interested in learning how the weights Wp 
are distributed among different eigenstates, how they de- 
pend on the initial temperature and quenching protocol, 
and how they compare to the weights of the eigenstates 
in the CE. Clearly, if Wp values approach the CE weights 
■^CE^"^'^ /^CE^ thermalization will take place. 

In Fig. [U we show the coarse-grained weights of the 
eigenstates of the final Hamiltonian in the DE [Figs.[lja), 
(He), and me)] and the CE [Figs, [^b), [^d), andUff)]. 
The results are obtained for quenches Aj = 4 to Ap — 0, 
and for three different initial temperatures. One can see 
there that, at the lowest temperature (T = 2), the distri- 
bution of weights in the DE [Fig. [TJa)] is very different 
from that in the CE [Fig. W[h)]. Wp exhibits a banded 
structure that is in stark contrast with the simple expo- 
nential decay seen in the CE. Remarkably, as the temper- 
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FIG. 1: (Color online) Density plot of the coarse-grained 
weights of energy eigenstates in the DE Wp (a, c, e) and 

in the CE Z^^.e'^^f' ^"^^^ (b, d, f). Results presented are for 
lattices with L — 24, A*' = 12 (half filling) and for quenches 
from A/ = 4 to Ap = 0. The values of the initial temperature 
are: T = 2 (a, b), T = 4 (c, d), and T = 8 (e, f). The color 
scale indicates the number of states per unit area. 



ature in the initial state increases, as shown in Figs. [Ijc) 
and md), the distance between the bands seen in the 
DE decreases and the slopes of the bands approach that 
of the CE. For even higher initial temperatures, higher 
than the gaps in the initial state (A cx Aj) [Figs.flje) and 
[Ijf)], the bands in the DE merge and the weights very 
closely resemble the canonical weights, i.e., the expecta- 
tion value of observables after relaxation in such quenches 
will be nearly thermal. It is not difficult to understand 
this because, in the T — > oo limit, the initial thermal 
state will have a completely flat (featureless) distribution 
of conserved quantities. Such a distribution of conserved 
quantities coincides with the one at infinite temperature 
in thermal equilibrium after the quench. Hence, the ini- 
tial state will essentially provide an unbiased sampling of 
the eigenstates that make the main contribution to the 
CE. 

Results for the reverse quench to the one in Fig. [1] 
namely, a quench from Aj = to Ap = are depicted 
in Fig. [2] Due to the energy gaps present in the many- 
body spectrum of the final Hamiltonian (generated by 
the presence of the superlattice potential), the distribu- 
tion of weights in both the DE and the CE exhibit bands 
separated by gaps A (x Ap. Once again, within each 
band, the weights in the DE and the CE are very differ- 
ent at low temperatures (T = 2 in Fig. [2]). However, it 
is also apparent in Fig. [5] that these weights become sim- 
ilar to each other as the temperature of the initial state 
increases. 

To better quantify the contribution of each part of the 
energy spectrum to the different ensembles, we have ex- 




FIG. 2: (Color online) Same as Fig. [T] but for quenches from 
Ai = to Ap = 4. 
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tracted out of Figs.[T]and[2]the energy density p{E). This 
quantity is proportional to the sum of the weights in a 
given energy window times the number of states in that 
window. In our plots, the results for p{E) are properly 
normalized such that the integral of p{E) over the full 
spectrum is unity. Results for that quantity, and T = 2, 
are presented in Fig. [31 It is remarkable to see that, for 
quenches from A/ = 4 to Ap = [Fig. [HJa) and (b)], 
the energy density in the DE exhibits a Gaussian shape 
independently of the filling of the system. 

In Ref. results for quenches starting from the 

ground state showed that only for half-filled systems and 
large values of Aj does a Gaussian shape develop in p{E) 
in the DE. Here we find that, at finite and not too low 
initial temperatures, the energy densities as well as other 
observables are qualitatively similar for different filling 
factors. This is in strong contrast with what happens for 
quenches starting from the ground state. Such a contrast 
is expected, as the ground states for different fillings may 
be qualitatively different, e.g., insulating at half-filling 
and superfluid away from half-filling, but those differ- 
ences are washed out with increasing temperature. In 
Fig. [21 note that the Gaussian like shape observed in the 
DE is slightly wider that in the CE. However, the width 
of both Gaussians is expected to vanish in the thermo- 
dynamic limit. 



A,=4 to A =0 




A,=0 to A =4 



(b) 



L=24, N=12 



. M iS JJ «■ w ^1 



0.3 



0.2 



(d) 



L=28, N=7 



- ,- » ° • 

I I ■ giFii^ii I 



0.1 





0.3 



0.2 



9. 



0.1 



-40 -20 




E 



20 40 



FIG. 3: (Color online) Energy density p{E) for two quenches 
at finite T — 2. (a, c) A quench from ^4/ = 4 to Ap ~ 0; 
(b, d) a quench from ^7 = to Ap = 4. (a, b) A half-filled 
system (A^ = L/2); (c, d) a quarter-filled one (A = L/A). In 
each panel, we report the results for p{E) in the DE and the 
CE. In (a) and (c), regardless of the filling, it is apparent that 
the energy distributions in both ensembles have a Gaussian 
shape. Hence, we also report results for a Gaussian fit p{E) — 

(V2^<5£)-ie-(^-^^''/(2*^') to each curve. 



The results for quenches from Aj = to Ap = A (de- 
picted in the Figs. [3]) also develop a Gaussian shape in 
both ensembles. However, this is less evident because of 
the presence of gaps in the spectrum. For these quenches, 
we also find that the results for the various observables 
studied here are qualitatively similar for different fillings 
when the temperatures are not too low. Because of this, 
in the remainder of the paper we focus on the half-filled 
case. 



IV. ENTROPIES 

We showed in the previous section that, for initial 
states at finite temperature, the energy distribution after 
quenches within integrable systems takes a Gaussian-like 
form. This is interesting because, for quenches start- 
ing from pure states that are eigenstates of the initial 
Hamiltonian, Gaussian-like energy distributions are only 
generically observed in nonintegrable systems [47| . Since 
the energy density is calculated through a coarse graining 
of the weights, in the following, we calculate the entropy 
associated with the DE and compare it to the entropy in 
the GE and the GGE. This allows us to better quantify 
the number of states contributing to each ensemble and 
to assess whether thermalization can take place in the 
system. 

In Fig. [31 we show the entropy per site with increasing 
L for the DE, GE, and GGE, and for quenches start- 
ing from different initial temperatures. There are several 
important conclusions that can be reached from those 
results, (i) For all entropies and quenches, S/L satu- 



A,=0 to A =4 




FIG. 4: (Color online) Entropy per site as a function of system 
size for various ensembles and initial temperatures, (a, c, e) 
Quenches from ^4/ = 4 to Ap = 0; (b, d, f) Quenches from 
Ai — to Ap =4. The systems are at half-filling and the 
initial temperatures are T = 1 (a, b), T = 2 (c, d), and T = 3 
(e, f). 
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FIG. 5: (Color online) Difference between S'gge and Sge per 
site as a function of L for quenches 7^ — > j4f = (a, c) 
and for quenches ^/ = — >■ /If 7^ (b, d). (a, b) Results 
obtained for a fixed initial temperature T = 1 but different 
values of Aj (a) and Af (b). (c, d) Results obtained for fixed 
Ai = 2 (c) and Af = 2 {d) but different values of T. All 
systems are at half-filling. 



rates with increasing system size, as expected since S is 
an additive quantity, (ii) For all quenches at finite tem- 
perature, Sde is a fraction of Sqe and S'gge, and our 
finite-size scaling analysis indicates that this will be the 
case in the thermodynamic limit, i.e., exponentially fewer 
states are involved in the DE compared to the other two. 
This is consistent with the findings in Refs. [l^, [2l|, \^ 
for quenches starting from pure states, (iii) Also, for 
all quenches, ^ge and S'gge are very close to each other. 
They can actually be seen to approach each other further 
as the temperature in the initial state increases. 

The lack of agreement between the entropy of the DE 
and that of the GGE does not mean that the GGE will 
fail to describe observables after relaxation. As shown 
in Ref. [l^, the great majority of eigenstates in the DE 
and generalized ensembles have very similar expectation 
values of few- body observables, i.e., independently of the 
number of states in each ensemble, the expectation value 
of the observables will coincide. What is of more interest 
in the remainder of the paper is how the GGE compares 
with the GE. An agreement between the two means that 
the integrable system would actually thcrmalizc. 

In Figs. [HI a) andlSfb), we plot a finite-size scaling of 
(S'ge — Sgge)/^ for systems with the same initial tem- 
perature but quenched between different values of Aj and 
Ap- In all cases, the difference is found to saturate to 
a finite value with increasing system size. Consistent 
with the findings in Ref. [2l| for quenches starting from 
the ground state, we do find that as the value of Aj in- 
creases, for quenches Aj ^ Ap = [Fig. O^a)], or as 
the value oi Ap increases, for quenches Aj ^ ^ Ap ^ 
[Fig. El^b)], the difference between the two entropies dc- 




FIG. 6: (Color online) (Sqe — Sgge)/L as a function of Aj 
for quenches Ai 7^0— >Af = (open symbols) and Af 
for A J = ^ Af t^O (filled symbols), with the initial 
temperature fixed. For each quench, results for two tem- 
peratures are presented. The dashed lines signal the power- 
law decay observed for large values of Ai{Af)- AJ^ for 
the quench Aj ^ ^ Af =0, and A'^'^ for the quench 
A/ = — > j4f 7^ 0. The systems are at half-filling with 
L = 32 in all cases. 



creases. 

In Figs. [5jc) and[5jd), we plot a finite-size scaling of 
(S'ge — S'gge )/i foi' systems with fixed values of Aj and 
Ap but for different initial temperatures. Similar to the 
results in Figs.lSfa) andlSjb), we find that the difference 
saturates to a finite value with increasing system size. 
However, that value decreases as the temperature in the 
initial state increases, and can become negligibly small. 

From the results depicted in Fig. [S] we arrive at the 
conclusion that, even though the difference between the 
entropy of the GGE and that of the GE becomes negligi- 
bly small if some control parameters (the values of Aj or 
Ap, and of T) are changed, it never vanishes in the ther- 
modynamic limit as long as those control parameters are 
kept finite and fixed. This means that the GE and the 
GGE do not become equivalent in the thermodynamic 
limit, and thermalization will not generally occur at a 
finite temperature. 

In Fig. [HI we show how the difference between the en- 
tropy per site in the GE and the GGE changes with in- 
creasing Aj ot Ap while the initial temperature is kept 
fixed. For both quenches, A/ 7^ — >■ Ap = and 
A/ = — > Ap 7^ 0, (Sge — S'gge )/-^ exhibits a power-law 
decay in the regime of large Aj{Ap). When the systems 
are quenched by switching off the superlatticc potential, 
the results for different T values are on top of each other 
once Aj becomes sufficiently large, and they decay at 
~ 1/A^. The exponent of this power law is the same 
found in quenches from the ground in Ref. (2]| . This 
is because the half-filled system in the presence of a su- 
perlatticc exhibits a gap A = 2Ai between the ground 
state and the first excited state, so as long as the tem- 
perature is much lower than the gap, the initial system is 
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FIG. 7: (Color online) (5ge — Sgge)/L as a function of T 
for quenches Ai ^ Af = (open symbols) and Ai = 
— > 7^ (filled symbols), for fixed values of Aj or Ap- 
For each kind of quench, results for three values of nonzero 
Ai{Af) are presented. Dashed lines are power-law fits to 
quenches Ai ^ 2 ^ Af = [(Sge - Sgge)/L ~ l/T®-"] and 
Ai = ^ Af = 2 [{Sge - Sgge)/L ~ in the region 

T > 10. The systems are at half filling with L = 32 in all 
cases. 



essentially in its ground state. (The ground state in the 
limit A/ — > oo is a trivial Fock state, namely, the prod- 
uct of empty and occupied single-site states.) For the 
quench Aj = ^ Ap 0, one can see in Fig. IH] that the 
results for different temperatures do not coincide with 
each other. However, they do exhibit the same power- 
law decay, (^qe — 'S'gge)/^ ^ ^Mfi independently of 
the temperature of the initial state. The vanishing of the 
difference between the entropy of the GGE and the GE 
as Aj{Ap) increases indicates that the ensembles become 
equivalent to each other and observables will exhibit a 
thermal-like behavior after relaxation. 

In Fig. [3 we study how (5ge — Sgge)/L behaves as a 
function of the temperature of the initial state for differ- 
ent quenches Aj ^ ^ Ap = and Aj = ^ Ap ^ 0. 
There one can see that the difference between the en- 
tropy per site in the GE and that in the GGE decreases 
with increasing T. The decay is power-law-like for high 
values of T, and is faster for quenches Aj ^ ^ Ap = 
than for quenches Aj = Q Ap 0. Wc have fitted 
our results for T > 10 to a power law in the quenches 
Aj = 2^ Ap = obtaining (5ge - Sgge)/L - l/T^-" 
and in the quenches Aj = Q ^ Ap = 2 obtaining 
(5ge - Sgge)/L ~ l/ri-99. 



V. CONSERVED QUANTITIES 

In this section, we study the expectation values of the 
conserved quantities /„ in the GGE and in the GE as one 
changes system parameters, similarly to what was done 
in Scc. lIVI for the entropy. The conserved quantities {/„} 
considered here are the set of L projection operators to 



the single-particle eigenstates of the final noninteracting 
fermionic Hamiltonian to which hard-core bosons can be 
mapped. By construction, the expectation values of the 
conserved quantities in the GGE are identical to those 
in the initial state, i.e., (/,i)gge = Tr[/9//„]. In the GE, 
they can be computed straightforwardly because the oc- 
cupations of the single-particle energy levels follow the 
Fermi-Dirac distribution (/„)ge = l/[l + e^'""'''/'^^^] , 
where e„ are the single-particle eigenenergies. 

In Fig. m we plot the expectation values of the con- 
served quantities in the GGE and the GE for differ- 
ent combinations of Aj{Ap) and T, for both types of 
quenches studied in the previous sections. Figures [DJa) 
and[51Jb) depict results for quenches with the same initial 
temperature but different values of Ai{Ap). When the 
half-filled system is quenched from Aj to Ap = 0, 
(-^n)GGE and (/n)GE are almost indistinguishable from 
each other and they are smooth functions of n [Fig.|SI^a)]. 
The picture for quenches Aj = Ap is very differ- 
ent [Fig. [Sl^b)] . In this case, the expectation values of the 
conserved quantities in the GE exhibits a discontinuity at 
n/L = 0.5, which is the result of the gap opened by the 
superlatticc potential. On the other hand, (/„)gge is a 
smooth function of n. However, the presence of the band 
gap does have an effect on (/n)GGE, as it tends to flatten 
out its values in the vicinity of the gap position (an effect 
that becomes more evident as Ap increases). This effect 
is only seen at finite temperatures, as in the ground state 
(In) GGE is the same between quenches from Aj and 
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FIG. 8: (Color online) Expectation value of the conserved 
quantities J„ = {In) corresponding to the nth lowest energy 
eigenstate in the single-particle spectrum. Results are pre- 
sented for quenches ^4/ 7^ — > Ap ~ (a, c), and for quenches 
^/ = — > 7^ (b, d). (a, b) Results for systems with the 
same initial temperature T — 1 but different values of Ai (a) 
and Af (b); (c, d) results for systems with the same Ai — 2 
(c) and Af = 2 (d) but different initial temperatures. In 
all panels, open symbols depict GGE results, and solid lines 
depict GE results. All systems are at half-filling with L = 36. 
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FIG. 9: (Color online) Integrated difference AJ as a function 
of L for quenches yl/ 7^ — = (a, c) and for quenches 
A/ = — tIf 7^ (b, d). Parameters are the same as in Fig. 
[5] All systems are at half-filling. 
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(Color online) A7 as a function of Ai, for quenches 
Af = G (open symbols), and of j4f, for quenches 
— >■ Af 7^ (filled symbols), for two values of T. 



FIG. 10; 

Ai^O^ 
Ai = Q 

Dashed lines depict power-law fits to the large Ai (A7 
l/Af) and Ap (A7 ~ 1/^/) results. Systems are at half- 
filling with L = 32. 



quenches io Ap when Aj in the former equals in 
the latter [2l|. Note that, for both kinds of quenches, the 
expectation values of all conserved quantities approach a 
constant (1/2) value as Ai or Ap increases. 

In Fig. [U we show the results for quenches A/ = 2 — 
Ap =Q [Fig. ^c)] 'atlA Ai = Q ^ Ap ^2 [Fig. ^d)] in 
systems with different values of the initial temperature. 
The overall picture is similar to that in the top panels 
when Ai{Ap) is changed. As the initial temperature in- 
creases, the expectation values of all conserved quantities 
approach a constant (1/2) value, while the specific fea- 
tures of each kind of quench are still visible. Namely, for 
quenches Ai = 2 ^ Ap = 0, the GGE and thermal pre- 
dictions are very close to each other and, for quenches 
Ai = ^ Ap = 2, a discontinuity is always seen in 
the GE prediction. Such a discontinuity in the latter 
quenches is absent in the distribution of conserved quan- 
tities in the initial state. 

In order to be more quantitative, and for comparison 
with the results obtained for the entropy, we calculate 
the relative integrated difference between the GGE and 
the GE predictions for the conserved quantities 

Scaling results for A/ vs L are presented in Fig. O They 
are qualitatively similar to those for the entropy differ- 
ences in Fig. [5l In all quenches studied, the differences 
between the GGE and the GE predictions for the en- 
tropies and the conserved quantities are seen to satu- 
rate to a finite value with increasing system size. Also, 
the differences between the GGE and the GE predictions 
for the conserved quantities decrease as Ai increases (for 
quenches Ai ^ Q ^ Ap = Q) and as Ap increases (for 



quenches Ai = Q ^ Ap ^ Q), as well as when T in- 
creases. 

The dependence of A/ on Ai{Ap), for a fixed system 
size and for two different temperatures, is depicted in 
Fig. [ini For both kinds of quenches, one can see that A/ 
vanishes as a power law in the regime of large Ai[Ap): 
for the quench from Ai ^ io Ap ~ A J decreases as 
A/ ~ 1/^/, while for the opposite quench it decreases 
as A/ ~' 1/ Ap. One should note that, as expected from 
the discussion in Sec. IIVI the behavior seen for ^/ 3> T 
is independent of T and identical to that in the ground 
state. The exponents of the power laws are also indepen- 
dent of T for both quenches and are found to be one-half 
those for (S'qe - S'gge)/^ vs Ai{Ap). 

A/ as a function of T is shown in Fig. [Til for three 
quenches Ai ^ Q ^ Ap = and for three quenches 
Ai = Q ^ Ap ^ 0. The behavior of A/ is once again 
similar to that observed for (iSge — Sqqe)/L vs T in 
the same quenches; A/ vanishes as a power law at very 
high temperatures. Power-law fits in the region T > 10 
yield AI - i/T^-O-i for quenches A/ = 2 -> = and 
AI - i/T^-oo foj. quenches Ai ^ ^ Ap ^ 2. Once 
again, the exponents are found to be one-half those for 
{Sge - Sgge)/L vs T in Sec. |TVl 

The preceding analysis clearly shows that only when 
the difference between the conserved quantities in the ini- 
tial state (in the GGE) and those in the thermal ensem- 
bles becomes negligible do the entropies and the ensem- 
bles themselves become equivalent. An understanding of 
this, as well as of the relation between the exponents of 
the power-law decays seen for (S'qe — Sgge)/L and AI, 
can be gained if one realizes that the entropy in the GGE 
and the GE can be written in terms of the occupations 
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T 



FIG. 11: (Color online) A/ as a lunction of T for quenches 
Ai ^ Q ^ Af = Q (open symbols) and quenches A/ = — >■ 
Af 7^ (filled symbols). Dashed lines are power-law fits in 
the region T > 10. We obtain A/ ~ i/r^ o* j-^^ quenches 
4/ = 2 ^ 4f = and A/ ~ 1/7100 for quenches Ai = 
— >■ Af = 2, respectively. All systems are at half-filling with 
L = 32. 

of the single-particle eigenstates as 

L 

5 = - ^[(1 - /„) ln(l - /„) + /„ ln/„], (17) 

where /„ = (/,i)gge for the GGE and /„ = (/n )GE for the 
GE. Once the occupations /„ in both ensembles are very 
close to each other, they are also very close to 1/2 (see 
Fig.©. One can write (/„)ge-1/2 = and (/„)gge- 
1/2 = e^*^^ (with \en\ <^ 1); it is then straightforward 
to show that 

L 

ri=l 

+ ^[(^^^^)^ - + 0[(e««^)« - (s«^)^]}, (18) 

which provides an excellent description of the results for 
{Sge-~Sgge)/L <^ 1 in Sec.lIVIand decays with a power- 
law exponent that is twice that for AI. 

VI. CONCLUSIONS 

In this paper, we have studied properties of ID lattice 
hard-core bosons (XY chain) after quenches that start 



from initial thermal states. In all cases considered, the 
quenches were generated by sudden changes in a super- 
lattice potential (a local space-dependent magnetic field 
in the spin language). We have shown that, in these in- 
tegrable systems after a quench starting from a finite- 
temperature state, the coarse-grained energy densities 
exhibit a Gaussian shape. However, the distributions of 
the weights are still qualitatively different between the 
quenched state and thermal states. This dissimilarity 
leads to an extensive difference among the entropies in 
the isolated system after the quench (the diagonal en- 
tropy), the GGE, and the thermal ensembles. On the 
other hand, as one (or both) of the two control param- 
eters explored here (the initial temperature T and the 
strength of the superlattice Aj or Ap) are tuned to in- 
finity, all ensembles become equivalent. The approach 
between the entropy of the GE and that of the GGE un- 
der such tuning was shown to be power law in the control 
parameter, independent of the quench protocol selected 
and of the values of the other parameters that were kept 
fixed. It is important to emphasize that, when all param- 
eters in the quenches were kept fixed and finite and the 
system size was extrapolated to infinity, the differences 
between the GGE and the GE results were seen to sat- 
urate at a finite value, i.e., these two ensembles do not 
become equivalent. 

We have also shown that such differences have their 
origin in the disagreement between the conserved quan- 
tities after the quench, which are determined by the ini- 
tial state, and the distribution of conserved quantities 
in thermal ensembles. By tuning the control parame- 
ters mentioned above to infinity, we have seen that the 
distribution of the conserved quantities in the quenched 
state approaches that in thermal equilibrium, which ex- 
plains why the generalized and thermal ensembles ap- 
proach each other under those conditions. However, we 
should stress that, in our particular systems of inter- 
est and quench protocols followed, thermalization only 
occurs when the control parameters are tuned to lead 
to completely flat distributions of conserved quantities, 
which is what happens when the initial temperature is 
infinite or when Ai{Ap) — > 00. 
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